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Abstract. 

We show how to create long range interactions between alkali-atoms in different 
|' hyper-fine ground states, with the goal of coherent quantum transport. The scheme 

uses off resonant dressing with atomic Rydberg states. We demonstrate coherent 
migration of electronic excitation through dressed dipole-dipole interaction by full 
solutions of models with four essential states per atom and give the structure of 
the spectrum of dressed states for a dimer. In addition we present an effective 
^ , (perturbative) Hamiltonian for the ground-state manifold and show that it correctly 

• ' describes the full multi-state dynamics. We discuss excitation transport in detail for 

a chain of five atoms. In the presented scheme, the actual population in the Rydberg 
& state is kept small. Dressing offers large advantages over the direct use of Rydberg 

i levels: It reduces ionisation probabilities and provides an additional tuning parameter 

for life-times and interaction-strengths. 
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1. Introduction 

Due to their remarkable properties PQ, Rydberg atoms emerge as a versatile tool in ultra- 
cold atomic physics. They can be used for diverse topics, such as quantum information 
[21 [3], atomic aggregates [4], the study of conical intersections [5], digital quantum 
simulations [6], electro- magnetically induced transparency [7], coherent population 
trapping [8], giant Kerr non-linearities [9] and ultra-long range molecules [TO] ITl] IT2] . 
since their long range interactions can be tuned over many orders of magnitude. However 
excited atoms in Rydberg orbitals are more vulnerable to ionisation [13] or spontaneous 
decay and more difficult to trap [HI [15] . 

Recently, several groups proposed to combine the advantages of both, ground-state 
and Rydberg atoms, through "Rydberg-dressing" schemes [161 El HI HH1 120] • The 
core idea is to use off-resonant laser coupling between ground- and Rydberg states to 
create eigenstates of the laser-coupled system in which a small Rydberg component is 
admixed to the ground-state. These resulting dressed states inherit some of the extreme 
properties of the Rydberg states, while preserving favourable properties of the ground 
state. Moreover, the additional degrees of freedom provided by laser Rabi-frequency 
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and detuning increase flexibility and promise potential shaping as well as dynamical 
control. 

Dressing has been extensively studied for van der Waals (vdW) interactions, 
generated when virtual transitions between levels of a two-atom system yield a distance 
dependent energy shift. Several articles have proposed Rydberg dressing techniques 
to induce this kind of interactions for ground-state atoms [HI [T71 [18], HH [201 El]. In 
this context it was shown that dressing ground-state atoms with Rydberg levels allows 
the large Rydberg vdW interaction energies to be reduced until comparable to typical 
energy scales of cold-atom traps, while simultaneously increasing ensemble life-times to 
those required for thermal equilibration effects [TBI [T7] . These are helpful features e.g., 
for realising a supersolid phase in Rydberg-dressed Bose-Einstein condensates [TBI [22] . 

Here, we demonstrate that the ideas behind dressing Rydberg vdW interactions can 
also be applied to resonant dipole-dipole interactions of neutral atoms. They occur when 
transitions between different states of the two- atom system are energetically resonant. In 
this case, the existence of an interaction potential is linked to electronic state transfer 
between atoms [23]. Dressing dipole-dipole interactions requires the use of two laser 
couplings and four electronic states (two ground- and two Rydberg-states), in contrast 
to one coupling and two states that suffice for vdW dressing. Hence, the scheme is more 
involved but also more flexible, offering a larger number of dressed states. It ultimately 
allows coherent quantum state transport between atoms in different long-lived ground- 
states over distances of 5-15 fim and for durations of many milli-seconds. 

Similar ideas have previously been applied to trapped Rydberg ions [M] . However, 
we are motivated by the possibility to generate potentials in an atomic many-body 
system that induce different kinds of (coupled) electronic and atomic motion dependent 
on the overall quantum state of the system. Useful applications of such potentials 
can be found in [251 E]- In many systems this goal is easier to achieve with neutral 
atoms than with ions, since the motion of the latter is typically dominated by their 
Coulomb repulsion or the indispensable trap. Furthermore the mapping of state-transfer 
interactions to normally non-interacting ground states by the laser dressing can be 
presented more clearly in the atomic case, where complications by trapping, mixing of 
internal and external dynamics and additional laser induced interactions can be ignored 
at least for free atoms. Finally, see Ref. [26] for dressing techniques targeting the control 
of atomic motion without linked excitation transport. 

Our results provide an additional handle on the time-scale of excitation transport 
and life-times, enable time-dependent control over hopping-strengths and can be used 
to vary the order of magnitude of dipole-dipole forces. Within some constraints [19J, 
the dressing also enables the use of well established trapping methods. As recently 
shown, dressed dipole-dipole interactions are an important tool for the realisation of 
atomic ring trimers [5]. In particular, the confinement of long-range interacting atoms 
on a ring, as required in Ref. [5], is greatly facilitated by the techniques presented here. 
Atomic ring trimers allow the detailed study of wave-function dynamics near and across 
conical-intersections. Dressed dipole-dipole interactions may also prove useful for the 
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study of exciton dynamics in Rydberg chains [U |27J . 

To demonstrate dressed dipole-dipole interactions, we employ numerical simulations 
of a model with four essential states per atom. In this case there are several dressed 
state manifolds, distinguished by the number of excited Rydberg atoms which they 
contain in the limit of vanishing dressing. We apply Van-Vleck perturbation theory to 
obtain analytical expressions for the induced effective interaction and to determine the 
parameter range within which the scheme can function. 

The paper is organised as follows: In section |2] we present the atomic states 
considered, describe the resulting Hamiltonian and outline the basic principles of the 
dressing scheme. In section [3] we consider an atomic dimer in detail, providing analytical 
expressions for the effective dressed interactions, diagonalising the full Hamiltonian and 
discussing different available dressing manifolds. In section H] we present two exemplary 
applications of our results, dressed ring trimers and exciton migration. We conclude 
in section [5] and discuss technical details of our calculations and the underlying atomic 
physics in the subsequent appendices. 

2. Dressed Rydberg aggregates 

2.1. Model 




Figure 1. (a) Schematic level diagram for laser dressing of a pair of ground state 
atoms with Rydberg states. The ground states \g), \ h) do not participate in inter- 
atomic interactions nor spontaneously decay on time-scales of interest. States | s), \p) 
arc highly excited Rydberg states, participating in binary long range interactions as 
explained in the text. Ground and excited states are coupled in a far detuned fashion 
as indicated in the diagram. The symbol 7 indicates the relevance of spontaneous 
decay, which we discuss in section T3. 51 (b) Implementation of the scheme sketched in 
(a) for 7 Li, using the indicated states to realise \g), \h), \s), \p). Also shown are the 
states energetically closest to | s ) , | p ) . For the hyperfine-split ground state, F denotes 
the total atomic angular momentum (nuclear, orbital and spin). See |Appendix B| for 
further details regarding the indicated transitions. 

We consider a collection of N alkali atoms, where the n'th atom is located at the 
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position R n . The distances between the atoms, Rki = |R/t — Rj|, shall be so large that 
interaction between two atoms can be neglected unless both of them are in a Rydberg 
state. In this case they experience strong dipole-dipole interactions. As discussed in 
[U |23j |27J [251 12H] the dipole-dipole interaction leads to exciton dynamics similar to that 
of molecular dye aggregates, molecular crystals or light harvesting systems [29J . 

For practical reasons we restrict ourselves to situations where the dynamics can 
be adequately described using only two "ground" states and two Rydberg states per 
atom. With the term "ground-states" we refer e.g., to two different hyperfme levels 
of the actual alkali-atom ground state. We shall denote the ground states by | g ) and 
| h ) and the Rydberg states by | s ) and | p ) . The latter correspond to the I = and 
/ = 1 angular momentum components of a Rydberg level with large principal quantum 
number v > 30. Schematically, the resonant dipole-dipole interaction U(TV) leads to a 
coupling between states | sp) and | ps ), whose matrix elements in the Hamiltonian take 
the form Udi P (R) \ sp){ps\. The transition strength is Ud ip (R) = —fi 2 /R 3 , where R is 
the inter-atomic distance and fi quantifies the transition dipole moment between the 
states | s ) and | p ) . We assume that all atoms are prepared in the m; = azimuthal 
quantum states and never acquire mi ^ 0. With atoms constrained in a two-dimensional 
(2D) plane, orthogonal to the quantisation axis, this ensures that there is no angular 
dependence of U (R) [231 12H] • Such a geometry covers both scenarios described later in 
this article. The interaction between atoms k and I reads 



To map the strong interactions Uki to the ground states, the atoms are irradiated with 
far-detuned dressing laser-fields that couple the ground-states and our two selected 
Rydberg states coherently. The relevant level diagram is sketched in Fig. [1] (a). In 
practice such a coupling is commonly achieved by two (or multi)-photon transitions 
[30] . so that the Rabi-frequencies VL SyP and detunings A S)P in the diagram have to be 
regarded as effective quantities. From the laser transition parameters, we assemble the 
effective dimensionless dressing parameters 



They are a measure of how "far-detuned" the laser coupling is, and will emerge as 
crucial quantities controlling the dressing, as it is the case for dressed vdW interactions 
[TEl HZ]. It is important that the state | g ) is directly coupled only to | s ) and | h ) only 
to | p ) . In order for this simple picture with four relevant states per atom to be valid, 
the detunings of both laser couplings have to be chosen such that all other transitions, 
coupling to further Rydberg states or ground state levels, are so far detuned that they 
can be safely neglected^! We show some realistic level diagrams for 7 Li in Fig. [1] (b) 

| If both couplings are realised by two-photon transitions, these considerations should include the 
virtual middle level. We require couplings | g ) O | m s ) O | s ) and | h ) -f-> | m p } •O- | p ) with uniquely 
assigned states | m s / p ). See |Appcndix B| for more details. 




(1) 



».p 



2A. 



(2) 
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to demonstrate how this constraint can be met in practice. Throughout the article, we 
will refer to the states | g ) , \s) as the "s-pair" and the states | h ) , | p ) as the "p-pair" . 

Using the four states introduced above as a basis for the single atom, an A-body 
basis state | k ) is written as 

\k) = \k 1 ...k N ) = \k 1 )®...®\k N ), (3) 

where kj G {g, h, s,p} describes the electronic state of the atom j. For example we write 
| ghs ) when the first atom is in state | g ), the second in | h ) and the third in | s ). After 
defining operators a^J, = \k n )(k' n \ with k, k! G {g, h, s,p} where n is the atom-index, 
the many-body Hamiltonian can be written as 

H = H + V, (4) 



with 



and 



h = -a s *i? - \ E & S + E Uni^tSL (5) 



n I 



n ^ ' n ^ ' 

We split if into two parts, Hq and V, to facilitate our later application of perturbation 
theory. The atom-light coupling in (j4j) has been treated in the dipole- and rotating wave 
approximations, as explained in Appendix A The exciton number operator 



^ = E(^ + 4f), (7) 

n 

measures the number of atoms in the p-pair of states. It is easy to see that commutes 



with H and hence the exciton number is conserved, see also Appendix C 

Of course equation (TjJ is a simplified description of a complicated multi-state 
system. For a more complete picture, one could include additional Rydberg levels 
adjacent to | s), \p), as shown in Fig.[T](b), van der Waals interactions of all four states, 
spontaneous decay and loss in the excited states. The latter arises as atoms in Rydberg 
states are lost due to auto-ionization or blackbody radiation [311 [32]. Importantly, 
parameter regimes can be found where these corrections are minor and (j4]) suffices, for 
example for the states and parameters shown in Fig. [TJ (b) and described in section 12.21 



2.2. Basic principles and exemplary parameters 

Consider an atomic sample such as a chain, prepared in states | 7r n ), in which the n'th 
atom is in the upper ground state | h ) and all others in the absolute ground-state | g ) 

kn) = \ggg-h-ggg). (8) 

Without dressing (Q s = Q p = 0) this is an eigenstate of the Hamiltonian (Tj0). If 
we ignore dipole-dipole interactions for the moment, an adiabatic ramp- up of the 
dressing strengths Q S)P will change this eigenstate into one with some admixture 
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of Rydberg population. Allowing then for dipole-dipole interactions finally causes 
excitation transport via processes of the schematic form: | gh ) — > \ sp) \ps) — > \ hg). 
Here — > symbolizes a transition of both atoms between the ground state manifold 
and the Rydberg state manifold due to the dressing lasers. Similarly ~> stands for a 
transition occurring due to resonant dipole-dipole interactions between Rydberg states. 
In section [3] we will confirm this expectation with numerical simulations and later justify 
and quantify it using perturbation theory. If the lasers are sufficiently off- resonant, the 
admixed fraction of excited state should remain small throughout. For this reason we 
can describe the process that we have just sketched as an effective interaction in the 
space spanned by the \n n ). We will refer to it as the ground-state manifold. For later 
reference we also define spaces spanned by 

\x n ) A = \ggg...p...ggg), | 7r n ) B = | hhh...s...hhh), (9) 

which we call gp (hs) single excitation manifold. 

Throughout this paper we focus on a specific realisation of our scheme, sketched in 
Fig. Q](b). We consider 7 Li atoms, whose 2 S\/2 ground-state has a hyperfine-splitting of 
Ahf = 800 MHz [33]. For the Rydberg state | s ) we pick a principal quantum number 
v = 80 and angular momentum I = state. For \p) we choose v = 80 and / = 1. The 
energy difference between the two Rydberg states is 4.5 GHz. We will assume detunings 
A s = A p = A = 50 MHz and dressing parameters a s = a p = a = 0.02 unless indicated 
otherwise. We will occasionally refer to a (without index), which implies 
For the above values of a and A, our Rabi-frequencies are Q s = Q p = 4 MHz. 

For the most part, we illustrate the dressing scheme through an atomic dimer with 
R = 6.5 jum inter-particle separation. For this distance, the strength of the bare dipole- 
dipole interaction is U = 92.5 MHz, just less than the doubly excited state detuning 
A s + A p = 100 MHz. The reason for this choice will become clear later. 

We stress that we have chosen these parameters to be specific, while our findings 
are more general. The constraints under which the parameters can be varied without 
invalidating the dressing scheme have been discussed in section [2~T1 All following results 
are expected to hold for any alkali species if these constraints are met. 

3. Dressed dimer 

For the principal demonstration of excitation transport through Rydberg dressing, we 
first treat the simplest possible case, the atomic dimer. 

3.1. Dressing dynamics in the ground state manifold 

In this section, we show dynamics from the initial state | ti\ ) = \hg). We consider two 
different ramps of the dressing couplings Q StP in the Hamiltonian: (i) a sudden jump from 
zero to their final value and, (ii) an adiabatic ramp over a small finite time T ramp = 0.1 
ms, effectively changing the initial state to 1 7f x ) = (| hg) + 0(a StP )[\ pg ) + | hs)])/J\f, 
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Figure 2. Transport dynamics of a Rydberg-dressed dimcr, (a) Populations n& = \ck\ 2 
of the two-atom states \gh) (violet) and \hg) (green) on a linear scale. The blue 
line is the total population in the ground state manifold. (b,c) On a logarithmic 
scale the differences of the two ramp-ups become visible: (b) sudden addition of 
the laser coupling: The singly excited states (red, yellow) are roughly suppressed by 
a 2 = 4 x 10~ 4 , indicated by the solid black line. Doubly excited states (blue) are even 
further suppressed, (c) adiabatic ramp-up of the laser coupling: singly excited state 
populations are precisely suppressed by a 2 . The expected spontaneous lifetime of this 
system is T c fj = (a 2 (r s _1 + t^ 1 )) -1 = 290ms as will be described in section [3751 



where M is a normalisation factor (see also section l3~2l) . |§|. 

We expand the time dependent quantum state as | \&(t) ) = Xlk Ck l k), where | k) 
are the iV-body basis vectors defined in ([3]). In this basis, we express the Hamiltonian in 
matrix form and solve the Schrodinger equation (SE) idt\^(t)) = H\ ) by matrix 
exponentiation. Using the four-state model introduced in section [27TI we employ the full 
Hamiltonian (jlj) without further approximations. 

For two atoms placed a distance R = 6.5 /jm apart, the resulting population- 
oscillations between the states | gh ) and | hg ) can be clearly seen in figure |2j On the 
logarithmic scale of panel (b), we see that also states containing one Rydberg excitation 
or more are populated, but their occupation is roughly suppressed by a 2 = 4 x 10~ 4 , 
indicated by the black solid line. The suppression of excited state population by a 
factor a 2 is well known from vdW dressing techniques [161 E32] an d a basic feature of 
off-resonant Rabi-coupling. As can be seen in figure [2} there will be many ground state 
population oscillations up to the expected spontaneous life-time of the dressed 2-atom 
state, which we estimate to about 0.29s as described in section [3751 Due to the dressing, 
the time-scale on which the interaction transports the excitation has been massively 
increased compared to the "bare" population oscillation period between Rydberg states 
\ sp) -H- | ps ), which would be T cx = 5ns. 

The comparison of panels (b) and (c) of figure [2] shows that the use of the more 
refined, adiabatically created, initial state has led to a more regular evolution of the 

§ To obtain dressed excitation transport from a state | tt\ ) , the ramp ought to be slow on the time-scale 
of the laser-coupling but fast on that of effective dipole-dipole transport. Thus Id rcss < Tramp < T ex 
with Xd r css = ^ 7T l^s.p = 0.25 fis and T ex ~ 4 ms (The value is read off figure [2] and defined as the time 
it takes for the excitation to oscillate from the first atom to the second and back.). 
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excited state populations, whose oscillations now precisely follow those of the two ground 
states. The population of states with a single Rydberg excitation never exceeds a 2 . 
Doubly excited states are further suppressed. However, this careful initial state creation 
is not actually required in order to observe dressed population oscillations, as we have 
just seen. 



3.2. Effective Hamiltonian 

The results of the preceding section already unambiguously demonstrate excitation 
transport through full numerical solutions of the problem. To further understand 
the population-oscillations more intuitively, we now consider the ground-state manifold 
spanned by 1 7Ti ) , 1 7r2 ) as the system of interest and its Hilbert-space complement as 
environment. The coupling between the two is treated as a perturbation. We use van 
Vleck perturbation theory [34], as outlined by Shavitt et al. [35] to derive an effective 
Hamiltonian in terms of the "system" only. This scheme conveniently takes care of 
degeneracies and generates an effective Hamiltonian well suited to describe excitation 
transport. A full analytical diagonalisation of (OS) is impractical even for N = 2. In 
contrast the perturbative results for N = 2 are intuitive and generalise straightforwardly 
to cases with more atoms. They are also valid over large ranges of parameter space as 
we shall show below. 

The aim of van Vleck perturbation theory is to find a basis that block-diagonalises 
the Hamiltonian (jlj) to a given order in the perturbation V. For V = this is achieved 
by the basis | k ) introduced in section [21 since Hq is already fully diagonal. The relative 
importance of V and Hq is governed by the dressing parameters a SjP . As the a SiP 
increase, the basis that block-diagonalises H becomes more complicated. Perturbation 
theory will thus be valid as long as the dressing parameters a StP are small. 

We will use the notation 

|k)^|k) = |k) + ^6 m |m), (10) 

m 

with b m ~ 0(a s>p ), i.e. |k) is a state whose leading component is |k) as long as a 
remains small and U is small (i.e. R is large). In this section we shall only be interested 
in the block of the effective Hamiltonian that governs the dynamics of the effective 
ground-state manifold | ffi ) = | hg ) and 1 7r 2 } — \gh). The details of the calculation are 
given in Appendix D here we merely present the results. In the basis \tti), \ n 2 ) one 
obtains 

H cS = (E 2 + E 4 )t + E(R), (11) 

where 1 is the 2x2 identity matrix. We have energy shifts (E2 + E4) and an explicitly 
distance dependent part E(R). The light shifts E 2 (-E4), corresponding to second 
(fourth) order in V, are given by 

E 2 = a 2 s A s + a 2 p A p , (12a) 
E 4 = -{atA s + a 4 p A p + a 2 s a 2 p {A s + A p )), (12b) 
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where we have used the dressing parameters defined in (J2J). These light shifts could 
allow optical trapping using the same fields that provide the dressing coupling. The 
position dependent part of the Hamiltonian, E(R), can be written as 

m = i f {R i i , da) 



with 



U 21 (R) W(R) 



W(R) = - (A s + A p ), (14a) 



where we have defined Cy = U^/(A a + A p ). As expected the states | tt\ ), 1 7r 2 ) acquire 
small excited state populations of the form: 

| 7Ti > = (\hg) + a p \pg) + a s \hs) + 0(al p ))/Af, (15) 

independent of R. Here M = ^Jl + a p + a 2 s is a normalisation factor. 

Let us briefly mention some special cases: For large inter-atomic separations, such 
that fy <C (A s + A p ), we can expand the expressions ( |T4j) and obtain W(R) ~ Woo = 
— a^«p(A s + A p ) and U(R) ~ c? s a^J\i. The shift Woo, which is independent of R, can 
then be merged with E±. For shorter distances, the factors 1/(1 — Uf 2 ) become relevant 
and diverge at [7 12 = 1, or equivalently at Uy(R) = (A s + A p ). We will see in the next 
section that this divergence can be traced to avoided crossings between the perturbed 
energy eigenvalues. We show in section 13.41 that ( 1141) fails only in a fairly narrow region 
around this avoided crossing. 

An interesting special case is A s = —A p . Here, the equations predict that 
the dressing effect on the ground-state manifold vanishes, which is confirmed by full 
simulations. We will however show in section 13.61 that this case has special appeal for 
dressing within the single excited state manifold. 

We also derive the contribution to if e ff of 6'th order in V in Appendix D.l We will 



make use of it in section 13.41 to illustrate the rate of improvement that can be achieved 
through higher-order terms of perturbation theory. 

3.3. Potential surfaces 

To get an overview of the consequences of dressing radiation beyond the ground state 
manifold, we now consider the full energy spectrum of the dimer as a function of 



inter-atomic separation. As we outline in |Appendix C[ the dimer Hamiltonian can 
be brought into block-diagonal form, with three blocks describing exciton numbers of 
N e = 0, 1, 2. Of these, only the single exciton block is non-trivial, hence we only show 



the eigenspectrum of its corresponding sub-matrix Mi (1C.2j) here. Explicitly we solve 



M 1 \m k ) = E k \m k ). (16) 
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Besides providing information on the possibilities of dressed excitation transport, the 
energy eigenvalues as a function of separation can also be viewed as adiabatic (Born- 
Oppenheimer) potential surfaces for motion of the atoms, as long as this motion is 
sufficiently slow [36] . 

The -R-dependent energy eigenvalues are shown in figure [3] for A s = 50 MHz and 
A p = 60 MHz, with other parameters as described in section 12.21 The slight offset 
between A s and A p makes it easier to grasp the structure of the resulting spectrum. At 
large distances R and for moderate dressing strength a, the eigenstates are essentially 
superpositions of two basis states, and have odd or even symmetry under exchange of 
atoms 1 and 2. These superpositions are listed in the caption of Fig. [3] (a), disregarding 
normalisation. Their energies are then mostly determined by the total detuning of 
the two superimposed basis states. As the distance R is reduced, the dipole-dipole 
interaction energy of the states consisting mainly of two Rydberg atoms (blue dotted- 
and red solid lines) becomes more prominent and eventually leads to avoided crossings 
between these and the other states of the spectrum. On the large energy scale of 
Fig. [3] (a), 0(A S)P ), the dipole-dipole shift of the | s~p ) ± | p~s ) states is most prominent. 
However also the other states with only a small doubly excited state component acquire 
a space dependent dipole-dipole potential, as can be seen in the close-ups, panels (b)- 
(d). The energy splitting of the states within the ground-state manifold, panel (b), is 
of order C(a 4 A SjP ) while that in the singly excited state manifolds, panels (c-d), is of 
order 0(a 2 A sp ). In panel (b) we additionally display the energy eigenvalues obtained 
from the effective Hamiltonian ( ITTj) . The shape of the potential is well reproduced. 

The dressed potentials shown in figure |3] (b)-(d) have at least two interesting 
features: (i) Their overall strength scales like ~ a 2 or ~ a 4 , it can thus easily be 
adjusted by choice of laser parameters and even manipulated time-dependently. (ii) 
The shape of the potentials can be modified beyond a simple attractive or repulsive 
form ~ ±i?~ 3 , due to the appearance of avoided crossings. For example the attractive 
branch of potentials in the ground-state manifold, shown in Fig. [3] (b), approaches a 
finite value for small separations R as can be seen from its continuation in Fig. [3] (a). 
This could avoid acceleration to high velocities during atomic collisions on the attractive 
potential, maybe reducing the probability of collisional ionisation [37j|27]. Light induced 
modifications of potentials that would be strongly attractive in the absence of dressing 
have already been discussed for the case with a single ground state in Ref. [26], and 
proposed as collisional shield for cold polar molecules in Ref. [38] . However, in our case, 
once the separation of the atoms is much smaller than the minimal value included in 
figure [HI which is R ~ 4/xm, we have to keep in mind that the simplified four-state 
model breaks down: The dipole-dipole shift of the dressed states then may become so 
large that they couple strongly to two-atom states that are not included in our model. 

Let us briefly discuss the case when the detunings have equal magnitudes. As A s 
approaches A p , the asymptotic \g~p) — \ p~g) energy curve (violet line) in Fig. [3] (a) 
is squeezed between the neighbouring ones until the states become degenerate. If we 
change the sign of one of the detunings and consider the special case A s = —A p , we 
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Figure 3. Energy spectrum of a Rydberg-dressed atomic dimer as a function of atomic 
separation R for A s = 50 MHz and A p = 60 MHz. We plot only energies in the single 
exciton sub-space, governed by the Hamiltonian (|C2|) . All states with odd particle- 
exchange parity correspond to solid lines, those with even parity to bullets, (a) View on 
the largest energy-scale set by the detuning. The two states with dominant components 
\ sp) ±\ps) for large R are most strongly affected by the dipole-dipole interaction. As 
a consequence, one of them undergoes avoided crossings with other states. Three 
avoided crossings can clearly be seen. The states are labelled in the legend according 
to their leading two-body content at large R. The state-content character changes each 
time a state undergoes an avoided crossing, so that for example the solid black curve 
corresponds to \gh) — \hg) at large R, but has become close to | sp ) — \ ps ) at small 
R. The vertical dashed line indicates the separation R = 6.5/im, chosen for dynamical 
examples throughout this article, (b) Zoom onto ground-state manifold. The crosses 
indicate the eigenvalues of the effective Hamiltonian l[TT|l . (c-d) Zoom onto singly 
excited state manifolds, the colour code is as in panel (a). Note the different energy 
scales of (b) compared to (c) and (d). 



obtain a particularly symmetric energy spectrum, shown in figure HJ As predicted by 
( fT4|) . the dressing-induced state transfer interaction vanishes between the states \gh) 
and \ hg). However symmetrical induced potentials can now be found in the singly- 
excited state manifold. The strength of these potentials scales as a 2 , different from 
those previously obtained in the ground-state manifold. We will discuss this further in 
section 13.61 
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Figure 4. The same as figure El but for A s = — A p = 50MHz. In contrast to 
figure [3] the ground-state manifold does not acquire a potential through the dressing. 
Symmetrical dressed potentials are obtained for the singly-excited state manifolds. 
The adiabatic potentials for the asymptotic \sp) ± \ps) states are cut-off at energy 
~ ±A due to avoided crossings. 



3.4- Performance of effective Hamiltonians 

In this section we revisit dressed dimer dynamics first presented in section 13.11 to apply 
the results of section [3721 Figure [5] shows a comparison of the evolution governed by the 
effective perturbative Hamiltonian ( ITTjl with that following from the full Hamiltonian 
01]). Already fourth order perturbation theory describes the effective induced hopping 
period fairly well, as can be seen in panel (b). When we move to the sixth order 
expressions, which can be found in Appendix D the agreement becomes even better: 
The deviation never exceeds 1% in the time-interval shown. In other simulations we 
found that perturbation theory agrees better with full calculations if a is small, as is 
expected. Even for parameters with poorer agreement between full- and perturbative 
dynamics (for example larger a or closer to an avoided crossing), we typically still 
obtain the primary features which make these induced potentials interesting: (i) We 
have persistent dressed hopping, (ii) the excited state occupancy is suppressed. 

Since we have included the dipole-dipole interaction in the unperturbed part of the 
Hamiltonian (jSJ), we do not a priori require it to be small compared to the detuning for 
the perturbation theory in section 13.21 to be valid. In practice we find that expression 
(II ip works well even for ~ A StP , as long as one stays clear of the avoided crossings 
at Uij = (A s + A p ). The large splitting of the eigenstates in the ground state manifold 
in the vicinity of the avoided crossing, see Fig. [3] (b), results in comparably fast 
population-oscillations. The regime where we can find dressed excitation transport 
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Figure 5. Comparison between exact dynamics using the Hamiltonian @ (a) and 
the errors resulting from using its perturbative approximations to fourth order (b, see 
and sixth order (c, sec (|D.4[0 for the transport dynamics of Rydbcrg-dressed 
dimcr from a dressed localized state. Panels (b) and (c) show the difference between 
the respective perturbation orders and the exact solution. Simulations begin from the 
state \hg) at t — — T ramp . Then, over a duration T ramp = lOO/xs, we linearly ramp up 
the Rabi-frcquencies f2 SjP to their final values. This approximately creates the state 
| hg ) at t = 0. The colour scheme is as in figure[2j The solid blue line in panel (a) shows 
the function p(t) = N(t)/N(0) = exp (— t / V e ff) , with r c g = 290 ms for orientation, see 
discussion in the next section. 



on time scales shorter than the life-time of the multi atom state hence is typically 
entered when Uij approaches A s + A p , but has not quite reached it yet. The effective 
population-oscillation-period in figure [5] is T cx = a _4 (l — f7^ 2 )/(2L r i 2 ) = 4.5ms. Here 
a~ A = 6.25 x 10 6 , and the correction factor is (1 — Uf 2 ) = 0.13 due to the vicinity of an 
avoided crossing. 

Let us reconsider the doubly excited state populations in figure [2j For atomic 
distances where the system is far from all avoided crossings shown in section 13.31 the 
doubly excited states would be suppressed by a 4 . This gets modified near these avoided 
crossings. For the distance R = 6.5 /xm as, the system is fairly close to an avoided 
crossing, thus the doubly excited states are more strongly populated. 

To conclude this section, we would like to remark on the position degree of freedom 
of the atoms carrying the interaction. In this article we treat it classically, assuming 
atoms with a precisely defined separation R. In practice, each atom will have a position 
uncertainty, for example due to the zero-point motion in a harmonic well. The resulting 
distribution of interatomic distances R can lead to a fast de-phasing of population 
oscillations such as shown in figure [5j For the example shown, if we would assume a 
separation uncertainty of ctr = R/20 the oscillations de-phase after 4 cycles. 

We however do not study de-phasing here, since it is not specific to the dressed 
interactions. If we considered direct dipole-dipole excitation transport, between the 
states | sp ) and | ps ) only, it would also de-phase after 4 cycles, which however would 
take place on a timescale of 10 ns due to the by a _4 (l — Uf 2 ) stronger interaction. Thus, 
whether or not de-phasing due to disorder poses a problem depends on the intended 
application of excitation transport and not on whether or not the interactions arise due 
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to dressing. 

In particular Ref. [5] discusses an exemplary application where it poses no problem, 
simulations presented in Ref. [5] fully include atomic positioning uncertainties. Also in 
references (27J ESI [28] we study systems with interesting combined excitation transport 
and atomic motion, with interactions that could be interpreted as arising via the 
techniques considered here. Again, de-phasing effects do not pose a problem. 



3.5. Spontaneous decay and loss 

The use of chains of Rydberg atoms for random walks and other aggregate physics 
[U |25] is limited by spontaneous decay and other incoherent loss processes from the 
Rydberg states. When using a dressing scheme, these decay processes are reduced by 
construction since the atomic population in Rydberg states is kept small. Throughout 
the paper, except section \3. 5.1\ we thus do not explicitly include loss in the presented 
dynamical calculations, but instead estimate the life-time of dressed Rydberg-aggregates 
based on their content of actual Rydberg population. For example the state | tt\ ) (|T5l) 
will be assigned a total effective decay-rate 7 c ff = a 2 s ^ s + a^7 P for two atoms. This 
scaling of 7 e g with the dressing parameters follows the same pattern as in vdW dressing 

PSE]. 

Most of our simulations treat excitation transport via v = 80 Rydberg states of 
litihum, which have life-times r s = 77 1 = 185.8 /is and t p = 7" 1 = 315.8 /is [31]. This 
yields r e ff = 7^ = 290ms for a = 0.02 and N = 2. The corresponding exponential 
population decay is schematically indicated in Fig. [5] (a). Note, that simple estimates 
of the life-time of a dressed system using the relations of this section and the results of 
section [3721 require excited state populations to be accurately described by perturbation 
theory. We have seen in figure [2] that this is the case, even fairly close to the avoided 
crossings in the spectrum. 



3.5.1. Density -matrix treatment To ascertain that the simple estimates of loss effects 
discussed above are valid for our cases, we also calculated results such as shown in 
figure [5] with the inclusion of loss and spontaneous decay. Consider the density matrix 

p = 5> k ,i|k)<l|, (17) 

k,l 

with basis-states as defined in fl3]). The prime on the summation symbol indicates that 
only states within the single exciton manifold (set bi in Appendix C) are considered. 
Let the time-evolution of p be given by the Master equation 

p = -i[H,p] + J2^ v iD k a \p] + TC[L,p], (18) 

ak 

where T> is the Lindblad superoperator V[a,p] = apa) — d)ap/2 — pa) a/2. For k = {s,p} 

define l{k) = {g, h}, then acts like | l{k) ){k \ in the electronic space of atom a and 
as unity elsewhere. These operators describe spontaneous decay. 
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The superoperator C has the form C[a, p] = dp — pa\ and the operator ( k \L\ k' ) = 
^Ryd(k)5kk'/2^- Here riR y d(k) counts the number of excited atoms in the state | k). The 
purpose of C is to describe loss from the Rydberg states to states external to our model 
[39], for example to further Rydberg states. For simplicity we assumed an identical loss 
rate V from | s ) and \p). 




t[msl t[ms] t[ms] 



Figure 6. The same scenario as figure El but calculated with a density- matrix 
formalism including spontaneous decay and loss, (a) All decay processes return 
population to the ground states \g,h), case (i) described in the text, (b) Some 
decay processes loose population from our system, case (ii) described in the text, 
(c) The same as (b), but on a logarithmic scale. The black line in (a,b) is the function 
p(t) = N(t)/N(0) = cxp (— t/r e fj), with r e g- = 290 ms, shown for orientation. 

We consider two cases: (i) All the population loss from the Rydberg levels described 
by the effective decay rates 7 SiP from [HI], including black-body effects, ends up in the 
respective ground-states. Here we set T = and 7s, P = 7 SjP . This is a worst case scenario 
for our purposes, as all incoherent processes within our four state model de-cohere 
excitation transport, (ii) We split r e ff from [HI] into the zero temperature component r 
and the black-body induced component r bbr via r^ 1 = Tq +t^ z . This yields r 0iS = 413.7 
/is, To iP = 1386 pis, Tbbr.s = 337.3 ps, Tbbr, P = 409.0 ps. We then assume that population 
returns to their respective ground-states with a rate 7s, P = (t 0| ( Sj p)) _1 , while black-body 
redistributed population is lost, either due to ionisation or because it leaves the single 
exciton manifolqjjj. We thus set T = r^ s . 

These calculations, shown in figure El verify that loss out of the system from the 
Rydberg states | s ) and | p ) has indeed no other consequences than an overall population 
decay. If spontaneous decay in the channels | s ) \g) and | p ) \h) is included, 
population oscillations additionally show a de-phasing. Neither effect is dramatic on the 
time-scales considered here. 

|| For example through a cascaded decay like | gp) —> ■ ■ ■ —> \ gg). 
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state T cff E/gg (3 = r cff /r cx 

\sss~.sp) a [/ 2 r C/£ = A,£ 

\ggg r --gp) i +(JV -i)^ « 2f/ A> jv-i^h> 

Table 1. Lifetimes and interaction strengths of various dressed iV-atom states: From 
top to bottom in the excited-state manifold, singly-excitcd-statc manifold and ground- 
state manifold. 

3.6. Dressing in singly excited state manifolds 

In section 13.31 we have pointed out that dressing induced dipole-dipole interactions can 
be created in the ground-state manifold (between the states \gh), \ hg)) or in two 
different excited state manifolds (between the states \sh), \hs) or \g~p), \pg))- The 
latter requires significantly different detunings A s and A p (such like A s = — A p ) to 
energetically separate the two different singly-excited state manifolds. We have explicitly 
verified population-oscillations for example between the states | hs ) and \sh). 

To know whether the ground or singly-excited state manifold would be more useful 
for a specific application it is important to consider the life-time of the dressed multi- 
atom state: The smaller a, the longer the life-time, but also the longer the effective 
hopping period T cx ~ 1/ (2U e g) required for a complete population-oscillation that 
transfers the excitation from one atom in the dimer to its counter-part and back. Here, 
the interaction U e g denotes the off-diagonal entry in the effective Hamiltonian of the 
corresponding manifold, e.g. U2i{R) in ( IT3l) for the ground-state manifold. 

Here we estimate the life-time of dressed states more roughly than described in 
section l3~5l assuming that ground-state atoms (\g), \h)) do not decay and excited 
state atoms (\s), \p)) decay with the same rate 7 = l/r . For a system of N atoms in 
various states we then again add the decay rates. 

For three different iV-atom states, the life-times determined in this way and the 
ratio of life-time and hopping period /3 = T c g/T cx are listed in table [TJ As long as 
a~ 2 < N, the ratio /3 for the singly-excited state manifold can be larger than for the 
ground-state manifold. 

Whether dipole-dipole interactions induced through dressing are advantageous over 
direct dipole-dipole interactions among Rydberg states, strongly depends on the problem 
at hand, as we discuss now for three examples. 

(i) Exciton migration: For the migration of a single or multiple excitons on a rigid 
chain, as studied in [I], the ratio (3 introduced above gives a direct measure of how 
many sites the excitation can traverse within the life-time of the whole chain. Evidently 
this measure is never improved by the present dressing technique, whose advantage for 
this scenario thus only lies in reduced ionization probabilites and simpler atom trapping. 
In section 14.21 we show an exemplary case where the overall life time in the presence of 
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dressing is long enough to exploit these features. 

(ii) Adiabatic transport: There are dynamical scenarios whose time-scale of interest is 
not directly given by T ex , such as the adiabatic entanglement and excitation transport 
that we reported in [25]. It can be seen that the ratio of the time-scale of transport 
dynamics to spontaneous life-time of the chain can be improved by factors of about 2 
by working in the dressed singly-excited-state manifold. 

(Hi) Ring aggregates: In section H~T1 we consider systems of three dipole-dipole interacting 
atoms trapped in effectively one-dimension on a ring. For such a construction dressing 
can be highly beneficial. 

4. Applications 

After choosing the dimer as the simplest case to work out the basic details of dressed 
dipole-dipole interactions, we now briefly present two exemplary applications of our 
results to larger systems. First, in section |4~T| we consider a dressed ring-trimer (N = 3 
atoms confined on a ring), a system of interest due to conical intersections of the 
adiabatic energy surfaces [5]. Then, in section I4.2[ we examine a dressed exciton on 
an atomic chain with N = 5. 

4-1. Flexible Rydberg ring-trimer 

In Ref. [5] we show that circular Rydberg trimers, consisting of three Rydberg atoms 
trapped in one dimension on a ring, exhibit interesting quantum dynamics near conical 
intersections (CIs). For practical realisations, however, dipole-dipole forces are typically 
too large for ring confinement. In this case the reduction of interaction strengths, which 
limits the parameter range where dressing is beneficial for exciton migration, turns into 
a benefit. As dressing also increases the system's life-time, it greatly facilitates the 
practical creation of dipole-dipole ring trimers. 

Consider three atoms tightly confined in one dimension on a ring, as sketched in 
Fig. [7] (a). Achieving such confinement is the first point in this section that is simpler for 
dressed ground state atoms than bare Rydberg atoms. On the ring, the bare interaction 
between atoms Uki is fully determined by the two relative angles 6*12 and #23 shown 
in the sketch. Adiabatic energies, defined by (fT6l) . then form 2D surfaces -Ex- (#12, #23). 
In the absence of laser-couplings, the three surfaces spanned by \pss), \sps), \ ssp) 
exhibit a conical intersection [ID], where two of them become exactly degenerate, at 
#12 = 023 = 27r/3. In the presence of laser-couplings, this feature is imprinted also 
onto surfaces spanned predominantly by ground states \hgg), \ghg), \ggh). Conical 
intersections are of great interest in particular in chemical physics, as they strongly 
affect the outcome of photo-chemical reactions [ID]. In cold atomic gases of Rydberg 
dressed atoms, they result in interesting non-adiabatic and geometric phase effects [5]. 
For a detailed description of this system and the concept of conical intersections we refer 
to Ref. [5] and references therein. 
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Figure 7. Adiabatic dressing dynamics according to (@| for temporally varying bare 
interaction strength Uki(t). (a) Geometric arrangement of dressed Rydberg ring trimer 
as discussed in Ref. [5]. The interaction strength is varied according to two different 
(classical) trajectories of motion of the ring-trimer as shown in the inset. For the 
solid path, symmetry forces the eigenstates of (j4]) to remain constant in time despite 
changes in Uki(t), this corresponds to straight crossing of the conical intersection of 
the (dressed) energy surfaces. For the red dashed path, dressed eigenstates change in 
time and are being adiabatically followed, corresponding to motion on the same energy 
surface in Ref. [5j. (b) Bare interaction strengths Uki(t) for the two paths, using the 
same coding. For comparison — 2A is also shown, where A = 24MHz is the excited 
state detuning, (c) Logarithmic plot of populations = |cfc| 2 , for the solid path in 
the inset of (a), using the same assignment of colours to states excitation number as 
in figure [21 Grey indicates triply excited states, (d) Populations n& = \ck\ 2 for the 
dashed path in (b), using a normal scale, (e) Difference between the exact populations 
in (d) and the expected dynamics according to the fourth order effective perturbative 
Hamiltonian (|D.7I) . 

For the purpose of the present article, it is more important to point out that the 
techniques worked out here are highly beneficial for the practical realisation of this kind 
of ultra cold conical intersections. This conclusion is reached through an extensive but 
technical survey of parameter space. We defer full details to a specialised publication, 
and present here only the core points: 

In order to realise the scenario sketched in Fig. [7] (a), initially at least three 
conditions have to be met: 

A. The dipole-dipole interaction energy should not exceed the transverse oscillator 
spacing u± to ensure one- dimensional dynamics. Our interest is in dynamics near the 
conical intersection, where the atoms form an equilateral triangle and are separated by 
d = \/3R, where R is the ring radius. Comparing the (possibly dressed) interaction 
strength with the strength of the ring trap, oj±, we have /i 2 /<i 3 a 4 < u±, thus preventing 
either too small ring radii R or to large interaction strength /i and hence principal 
quantum numbers v. 

B. The motion should be adiabatic, except close to the conical intersection. We can 
estimate the order of magnitude of the CI transit time Tqi from the time-scale set by 
the classical equations of motion 9 = —(dV(6)/d6)/(2MR 2 ), which holds for a dimer of 
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angular separation 9. This time scale is Tqi = 2 3 / 4 V / Mi? 5 / 2 /(/i« 2 ). Adiabatic dynamics 
must fulfill Tci T cx , with T ex as in section 13.61 

C. The life-time of the dressed three-atom system should be longer than the duration of 
the motion, thus r c ff > Tci- 

Attempting to fulfill the inequalities (A-C) for various atomic species leads to two main 
conclusions: (i) Lighter atoms are favourable, (ii) Without dressing (a = 1), principal 
quantum numbers v would be too high to avoid strong effects of black-body radiation. 
If we consider a dressed system however, one more constraint comes into play: 

D. The bare dipole-dipole shift without dressing must be smaller than the doubly excited 
state detuning A s + A p , at the closest approach of the atoms. This is to stay clear of 
avoided crossings such as seen in figure [31 

A set of parameters where these conditions are met is: v = 100, R = 9.8 /im, 
A s = A p = 24 MHz and a s = a p = 0.150. Now we consider the dynamical evolution of 
an eigenstate of (J4]), when the interaction parameters U^i = U^iif) are varied in time in 
a manner representative for near CI dynamics. To this end we extracted the trajectory 
0(f) = (6) from full quantum mechanical solutions of the time-dependent Schrodinger 
equation in Ref. p)fn . shown in the inset of Fig. [7] (a) as dashed line. A second trajectory 
is taken slightly offset (dotted line). We find that the quantum state dynamics in both 
cases follows if e ff (see equation (1D.7|) ) well, justifying the simplified model used in [5]. 
Since the energy-spacing E p d (see figured]) for v = 100 is E p d = 298 MHz, we could have 
chosen a much larger detuning, yielding even better agreement between the perturbative 
and exact solutions. 

As we have seen in this section, dressing is threefold beneficial for atom trapping in 
the construction of conical intersections in dressed Rydberg ring trimers: (i) It greatly 
widens the available scope for quasi-lD trapping of the atoms, (ii) reduces the involved 
interaction strength such that trapping can be realistically considered in the first place 
and (iii) simultaneously extends the available life-times. 

4-2. Dressed excitation transport on long atomic chains 

We discuss results for a linear chain of N = 5 equidistant atoms in the following, 
considering the same physical situation as in section 13.41 hence the spacing of atoms in 
the chain is also AR = 6.5/im. We create the initial state adiabatically as in section l3~4l 
For this section only, we will define our excitation to be the state \g), such that 
e.g. 1 7Ti ) = | ghhhh). The corresponding dressed states have a longer life-time for larger 
chains. 

The dressing-induced excitation transport and a comparison of exact and 
perturbative evolution is shown in figure [HJ It can be seen that the excitation is 
transported over many sites of the chain on time-scales shorter than the expected 

These are the parameters employed in [5], in that article the detuning and Vb are (d*) were given as 
angular frequency A u — 2ir x A. 

+ These simulations only contained the three effective ground-state energy surfaces. 
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Figure 8. (a) Dressing-induced excitation transport in a chain of 5 atoms, according to 
the Hamiltonian (|D.7|) . Visualized is the population of states \TT n ). (b) Quantitative 
details of the same data as in (a). Sites 1 and 5 are shown in (magenta, green) 
respectively, the approximate population loss is shown as blue dashed line, (c) Absolute 
difference between the evolution according to the fourth order effective Hamiltonian 
(|D.7|) and the full 5-body Hamiltonian (gj. 

spontaneous life-time of the dressed 5-atom state (about r e ff = [a 2 ((N— l)/r p +l/r s ] _1 = 
0.14 s). Due to the dressing, the time-scale on which the interaction transports the 
excitation has been increased by a factor of a _4 (l — U12) = 8.4 x 10 5 , compared to the 
direct use of excited Rydberg states. 

In section 13.61 we argued that dressing typically does not favourably alter the 
possible number of excitation hops, given by the ratio (3 = r c ff/T ex of life-time and 
hopping period. Nonetheless the results of this section show that dressed excitation 
migration within the available life-time is possible. Despite the detrimental effect on f3, 
the dressing still enables time-dependent control of the transport, and, as in the previous 
section, simplifies trapping efforts. 

5. Conclusions and Outlook 

We have shown that the dressing of alkali atoms in the ground-state with Rydberg- 
excitations can yield dipole-dipole interactions with excitation-transfer for the ground- 
state atoms. This generalises existing results for dressing with Rydberg van der Waals 
interactions, that do not entail excitation transport. The scheme proposed here is an 
adaptation of a similar method, introduced in the context of trapped Rydberg ions [24] , 
to neutral atomic systems. It makes use of two effective laser couplings to Rydberg- 
states, instead of the one required for van der Waals dressing. 

We find that the dressed dipole-dipole transition matrix element between ground 
state atoms scales like a 4 with the dressing parameter a, while the life-time of the dressed 
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state increases like ct~ 2 . Dressed vdW interactions show the same scaling behaviour. 
We demonstrate in Ref. [5] that this re-arrangement of life-times and interaction time- 
scales allows the study of conical intersections using Rydberg-dressed atoms on a ring. 
Beyond this, the dressing scheme will enable a larger variety of trapping techniqes for 
the atoms and reduce the likelihood of ionisation. As further example of the flexibility 
afforded by the dressing, we show that it is possible to induce effective transport of a 
single Rydberg excitation within a chain of ground state atoms. The amplitude of this 
process scales like a 2 . 

Whether the effective interactions discussed in this article are more useful than the 
direct use of dipole-dipole interactions depends on details of the situation one studies. 
In Ref. [5] we provide an exemplary scenario where they can be highly beneficial. 
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Appendix A. Atom-light interactions 

In this section, we derive the rotating-wave Hamiltonian (Tj0) from more fundamental 
expressions for atom-laser interactions, reviewing text-book material, see e.g. [H]. The 
Hamiltonian for an atom in the presence of an electric field is 

H' = Wo-er -E(r 0j *). (A.l) 
This expression is valid in the dipole approximation k ■ r 1, when the wavelength of 
the incident light is much larger than typical atomic distances. The operator Wq is the 
full atomic Hamiltonian in the absence of any external field, e is the electron charge, 
r the position operator and E(ro,t) the electric field at the atomic nucleus, which is 
treated classically. 

Consider for the moment a single atom only, with the four relevant states | g), \h), 
\s), \p) and a level scheme as shown in Fig. [TJ (a). These states are eigenvectors of the 
free atomic Hamiltonian Wq, thus Wo\k) = HEf. with k G {g, h, s,p}. We add two 
coupling fields E = E s +E p , with E SiP = S SjP z cos {v s>p t) = £ StP (t)z. Here z is a unit vector 
in the ^-direction. For any two states | a ), \ b ) we define matrix elements £i a t, = e{a\z\b). 
A general quantum state in the chosen subspace is | ^(t) ) = J2 n e{ g h s p} c «l n )■> with 
equation of motion 
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* For each of the two transitions that we require, we only consider a single light-field instead of 
the multiple lasers that would typically be employed for a multi-photon transition. A more detailed 
treatment does not add qualitatively new features. 
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Here we have used atomic units. Now we change the variables c m of the coefficient vector 
to d m = c m e tuJmt , where u m are arbitrary constants. If the coefficient vector c evolved 



according to ij^c = Mc for some matrix M, d evolves according to ij^d 



M' = EME* 
and fi mm = ^m- 

Using the specific transformation vector u> 
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(A.3) 



where = E* 



Eg, Up 



2 fthp^pf 

E p — E h and f±(cu, v) stands for e ±i{uJ+u ^ + e ±«"-")t. Under the 
rotating-wave approximation all terms e ±t ^ s ' p+Us ' p ^ t are dropped, as they are oscillating 
too rapidly to have any effect. We also introduce the detunings A StP = v sp — u StP and 
Rabi-frequencies Q s = —Hg S £ s and Vt p = —Hhp£ P - 

Finally we perform one further variable-change, using uj = (0, A s , 0, A P ) T to reach 



(A.4) 



The Hamiltonian in (IA.4j) is the basic building block of (jl]), as long as only a single 
atom is concerned. Even for more than one atom and with the inclusion of dipole- 
dipole interactions, the above procedure can be followed. Instead of the single atom 
energies in the first transformation vector uj, we would employ many-atom energies. 
Since the dipole-dipole interaction only couples energetically degenerate states like | sp) 
to \ ps), all complex phase-factors can finally be eliminated to arrive at a many-body 
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Appendix B. Lithium energy levels 

In practice the scheme displayed in figure [1] is complicated by selection rules and 
transitions are constrained by the availability of laser sources. The coupling between the 
states | g } <H- | s ) can be realised with now broadly established two-photon excitation 
schemes. These transitions are typically near resonant with some auxiliary middle 
level, such as | 2p 3 / 2 ) in figure IBll in order to enhance transition amplitudes. For 
the right choice of parameters, coherent coupling between \g) = |2si/ 2 ,-F= 1) and 
\s) — |80si/ 2 ,-F = 1) is achieved, see e.g. [12]. For the coupling between the states 
| h) •H- \p) two-photon transitions are forbidden due to the selection rule Al = 1. 

This necessitates one- or three photon coupling schemes, examples for which are 
included in figure IBll Let us briefly discuss both options: 

(i) Single photon XJV -transition: Direct transitions from ground- to Rydberg levels [43] 
suffer from extremely small dipole-transition matrix elements. These scale like z/ _3//2 
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Figure Bl. Several alternative transition path- ways to realise the overall effective 
coupling between the ground- and Rydberg states. The | g ) o | s ) coupling proceeds 
via a standard two-photon transition, while the | h ) -f-> | p ) coupling requires either 
a UV-transition [43] or a three-photon transition [44] . both of which are more 
cumbersome. This diagram is intended as a rough guide only, hence fine- and hyperfine 
structure are omitted for all excited states. Intermediate state energies are taken from 
05], Rydberg energies calculated with simple quantum defect theory as in [4H] . 



with the principal quantum number pp. From the matrix element |( 13^3/2 \r\ 2s\/2 ) \ = 
0.036a |47][J we extrapolate M = | ( 80p 3/2 \r\ 2s 1/2 >| = 0.0024a , where a is the Bohr 
radius. Assuming P = 0.09 mW of laser-power, focussed to a waist w — 10 /xm, one 
achieves a Rabi-frequency Q = 4 MHz as employed in this article. The intensity at the 
focus is / = 2P/ (ttw 2 ), resulting in a Rabi-frequency f2 = AieS/H, where S = a/2// (ceo) 
is the electric field, and e, c as well as eo fundamental constants. Laser light at the UV 
wavelength is typically created through higher harmonic generation, reducing flexibility 
and making it more challenging to achieve the required power. 

(ii) Three photon transition: One would probably employ three different transitions, as 
indicated in figure IBll exploiting near resonant auxiliary levels to enhance the coupling 
amplitude. In principle it would also be possible to utilise three photons stemming 
all from the same laser. For Li such three-photon, single colour, schemes still can 
proceed near resonant with the intermediate | 2pi/ 2 ) state |H]. Consider for example 



(j We calculate matrix elements M. from oscillator strengths fik and transition energies AE using 
A4 = \Z3e 2 h 2 fik/ (2m e AE), m e is the electron mass. 
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the excitation chain sketched in yellow and dark blue in figure IB1[ and denote the Rabi- 
frequencies and detunings of the three transition as Q^, A& with k = 1, 2, 3 from lower 
to higher energies. Since both intermediate states decay on a timescale of 10 ns, they 
have to be far detuned with ai = Q\/(2A 1 ) and a 2 = r2 2jC ff/(2(A 1 + A 2 )) of the order 
of afc < 10~ 3 for k — 1, 2. Here f2 2i efr = a\Q 2 . This makes sure that the effective decay 
rate of the intermediate levels is at least not larger than that due to the Rydberg state. 
If we now consider the effective three-photon Rabi frequency: 

ficff = 4A l( A 1 + A 2 ) = (R1) 
we would require O3 >4x 10 3 GHz to reach Q c s = 4 MHz as employed here, rather 
hard to achieve on the weak Rydberg transition. There may though still be applications 
where much smaller effective Rabi-frequencies Q p are sufficient and the three-photon 
scheme has advantages over the UV transition. 



Appendix C. Multiple excitons 



The Hamiltonian given in (T4J) is the iV-body generalisation of the results of Appendix A 
with added dipole-dipole interactions. It conserves the exciton-number, represented by 
the operator ([7]). Consequently the Hamiltonian has a block-diagonal structure, with 
each block describing a given number of excitons. For N = 2 atoms we explicitly have 
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with 

\ s 9}; 

I hp), 



The bases 

ilgg), \gs) 

b 2 = {\hh}, 

this Hamiltonian in more detail. 



respect to which the three blocks are written are b = 

\ ss )}i b i = {l5 ,/i )> \9P), \hg), \hs), \sh),\sp), \pg), \ps)}, 
\ph), \hh)}. In section [3731 we study the nontrivial part, Mi, of 



Appendix D. Van-Vleck perturbation theory 

The basic goal of Van-Vleck perturbation theory was outlined in section 13.21 to find a 
basis that block-diagonalises the Hamiltonian (T4j) to a given order in the perturbation 
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V. This appendix supplies the details necessary to understand the origin of the results 
presented in section 13.21 and Appendix D.l We first partition our unperturbed basis 
| n ) into two sets V and Q, the first of which shall span the "system" space of interest 
and the second is its complement in the full Hilbertspace (the "environment"). The 
specific basis, | n), used for the definition of system and environment is given in fl3]). 
This basis is also an eigenbasis of H . Then we can construct projection operators on the 
system subspace P = Yl<pev I ^ )( ^ I ari( ^ ^ s complement Q — 1 — P. Here we consider 
the example where V = {\n n ),l < n < iV}, thus our system space is the ground-state 
manifold. We then have P = J2 n I ) ( I • 

With respect to these partitions the Hamiltonian matrix, or similarly any other 
operator, can be divided into four blocks PHP, PHQ, QHP, QHQ. Out of these, 
we assemble a block diagonal part Hp = PHP + QHQ and a block-off-diagonal part 
Hx = PHQ + QHP. Since Ho is diagonal in its eigenbasis, we have Ho = Hq + Vd, 
Hx = Vx- We now aim to find a unitary transformation T that yields a block diagonal 
Hamiltonian H = T~ 1 HT to a given order in V, thus Hx = QHP + PHQ = 0. The 
effective Hamiltonian in the space of interest that we seek, is then given by the block 
H ea = PHP. 

Next, we express the unitary transformation operator T as T = e G with G = —G\ 
We then introduce the condition Go = and Gx = G for the block diagonal and 
off-diagonal parts of G. Other choices are possible, the present one distinguishes 
the van Vleck procedure from other related schemes [35]. For all operators in the 
problem, we write a series expansion in orders of the perturbation V, most notably 

n = H + ZZi w {n \ T = ZZi T{n \ g = ZZi G{n) - 

As outlined in [35], it is then possible to obtain recursive relations for [H , £J 
One extracts the full matrix G^ n > from [H , G^ 1 '] as follows. Let us introduce the notation 
| i ) with latin indices to denote an eigenvector of H within Q and | a ) with greek 
indices for one within V. Since we demand Go = we only require matrix elements 
like ( i | G^ \a) or ( a \ G^ \ i ) in order to know the whole matrix form of G^ . Now we 
define the resolvent operator. 

f = E M > (D.l) 

which fulfills: 

O\a) = -RS ) [H ,O]\a). (D.2) 

for any operator O. Substituting O = G^ in flD.2j) and multiplying from the left with 
(i| we obtain (i\G^\a) and also (a\G^\i) = -(i\G^\a)*. 

From knowledge of G one can infer the expansion orders through commutation 
relations. A series expansion of T = e G to the required order finally yields the perturbed 

ff A crucial step is to rewrite H = e~ G He G as series of increasingly deeper nested commutators 
using the Baker-Campbell-Hausdorff formula. By segregating expansion orders in V, the equation is 
then turned into a recursive set of equations for the [Ho,G^]. An elegant method to perform this 
cumbersome routine is outlined in [35], to which we refer for further details. 
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(D.3) 



Wc refer to Rcf. [351 for all further technical details. 



Appendix D.l. Higher orders and atom numbers 

The effective Hamiltonian of the dimer in the ground-state manifold spanned by 1 7^2 ) 
is given in section 13.21 up to fourth order, with odd orders vanishing. Here we 
display the corresponding expression of sixth order, Hq\ where we use the definition 
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{uh - 1) 2 

Recall thatA = A s + A p and U 12 = U 12 /A. 

Finally, for both cases shown in section HI the effective Hamiltonian of section 13.21 
has to be adjusted for N > 2. In terms of the basis | ft n ) we obtain: 

H e s,ij = (E 2 + E 4 )5ij + Uij(R), 
= (N - l)a 2 A s + ajA,, 



E 2 



-[(N - l)atA s + a*A f 
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- A p ) 8ij 
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(D.8) 
(D.9) 
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